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ABSTRACT 

We study the Sunyaev- Zeldovich effect for clusters of galaxies. The rate equations for 
the CMB photon distribution function are studied in three Lorentz frames. We clarify 
the relations of the SZ effects among the different Lorentz frames. We also derive 
analytic expressions for the photon redistribution functions. These formulas will be 
extremely useful for the analysis of the forthcoming observations. 
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1 INTRODUCTION 



The Sunyaev-Zeldovich (SZ) effect l|Zeldovich and Sunvaevl 
Il969l : ISu'nvaev and Zeldovich! 1 19801 ). which arises from the 
Compton scattering of the cosmic microwave background 
(CMB) photons by hot electrons in clusters of galaxies 
(CG), provides a useful method f or studies of cosm ology. 
For the reviews , for e xample, see iBirkinshawl (j 19991 ) and 
ICarlstrom et all (|2002l ). The original SZ formula has been 
derived from the Kompaneets equation (|Kompanccts fl956l ) 
in the non-relativistic approximation. Ho wever, X-ray obser- 
vations, for example, lAllen et al.l l|2002l ) have revealed the 
existence of high-temperature CG such as fcsT e ~20keV. 
For such high-temperature CG, the relativistic corrections 
will become important. 

On the other hand, it has been known theoretically for 
some time that the relativistic corrections become significant 
at the short wave length region A < 1mm. In particular, the 
recent report on the first dete ction of the SZ effe ct at A < 
650/im by the Hershcel survey (|Zemcov et al.|[2010t ) seems to 
confirm the relativistic corrections. F urthermore, new gen- 



eratio n observations, for example, by IPlanck Collaboration! 
(|201lT ) are carrying out systematic studies of the precision 
measurements on the SZ effect. Therefore, reliable theoreti- 
cal studies on the relativistic SZ effect at short wave length 
will become extremely important for both existing and forth- 
coming observation projects. 

The theoretical studies on the rela tivistic correc- 
tions have bee n done by sev eral groups. IWrightJ (|1979| ) 
and Rephaeli (|Rephaelil Il995l ) hav e done the pioneering 



work to t he SZ effect for the CG. Challinor and Lasenbvl 
( 1998j) and lltoh. Kohvama and Nozawah 19981 ) have adopted 



a relativistically covariant formalism to describe the 
Compton scattering process and have obtained higher- 
order relativistic corrections to the thermal SZ ef- 
fect in the form of the Fokker-Planck expa nsion ap- 
proximation. iNozawa. Itoh and Kohvamal |l998) have ex- 
tended their method to the case where the CG is 
moving with a peculiar velocity /3 C with respect to 
the CMB and have obtained the relativistic correc- 
tions to the kinematical SZ effect. Their results were 
confirmed by Challinor and Lasenbv |l999T) and also by 
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ISazonov and Sunyaev ( 19981 ). Itoh. Nozawa and Kohvamal 
(2000) have also ap plied the covariant formalism to the po- 
larization SZ effect (|Sunvaev and Zeldovichlll98ll ). The im- 
portance of the relativistic corrections is also exemplified 
through the possibility of directly me asuring the clu ster tem- 
perature us ing purely the SZ ef fect |Hansenll2004l ). On the 
other hand. lChluba et all l|2012l ) also studied the relativistic 
corrections on the SZ effect based on an alternative deriva- 
tion of the Boltzmann equation. They reported the clean 
separation of the kinetic and scattering terms which usu- 
ally mixed in the CMB frame due to the retardation effect 

In the present paper, we give formal relations for the 
rate equations among different Lorentz frames, namely, for 
the CMB frame, CG frame and a general observer's (OBS) 
frame. Starting from the Lorentz invariant Boltzmann equa- 
tion for the photon-electron scattering in the CMB frame, 
we derive the relation formula between the CMB and CG 
frames. First, we derive the expression for the rate equa- 
tion of the photon distribution function in the CG frame. 
Then, we transform the rate equation to the CMB frame. We 
compare the obtained expre s sion w ith the previous work of 
INozawa. Itoh and Kohva ma (1998]). We will show that two 
calculations agree completely, if we multiply the Lorentz 
factor l/7c to the result of INozawa. Itoh and Kohvamal 
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ll 19981). The reason for the fac t or is because N£ was used 



Nozawa. Itoh and Kohvamal (1 19981 ) instead of N e for the 



electron number density in the CMB frame, where N e and 
Aff are the electron number densities in the CMB and CG 
frames, respective ly. The factor is con sistent with the com- 
ment indicated bv lChluba et al.l (|2012h . This clarifies the re- 
lation of the rate equations between the two Lorentz frames. 
We will also extend the rate equation in the CG frame 
to the OBS frame. Finally, we derive analytic expressions 
for the redistribution functions Pe, c (s c ,P) for I = 0,1,2. 
With these expressions, 5-dimensional integrals for solving 
the rate equations are reduced to 2-dimensional integrals, 
which will save the numerical computation time drastically. 
These analytic expressions will be extremely useful for the 
analysis of the forthcoming observations. 

The present paper is organized as follows. In Section 2, 
we derive the formal relations for the rate equations of the 
photon distribution functions among three Lorentz frames. 
In Section 3, we derive the rate equation in the CG frame. 
We also express the rate equation in the operator represen- 
tation and in the Fokker-Planck expansion approximation. 
In Section 4, we transform the rate equation in the CG frame 
to the CMB frame. We will compare the present result with 
the previous works. In Section 5, we will transform the rate 
equation in the CG frame to the OBS frame. In Section 6, we 
derive analytic expressions for the redistribution functions. 
Concluding remarks will be given in Section 7. Finally, we 
show the relations for the electron number densities among 
the three Lorentz frames in Appendix A. In Appendix B, we 
give the normalization of the photon redistribution function. 



2 GENERAL FORMALISM 

In the present section, we discuss the formal relations of the 
rate equations for the CMB photons among three different 
Lorentz frames, namely, the CMB frame, the CG frame and 
the OBS frame. 

Let us start with the following Lorentz invariant ex- 
pression of the Boltzmann equation for the electron-photon 
scattering in the CMB frame: 

dn(ui) r> „ , » 
uj „\ ' + k ■ Vn(cj) 



Ot 



d 3 p d 3 p' d 3 k' r4 



j 5 (p + k — p — k ) a X 



(2tt) 3 J E E' 2oj' 
[n( U ){l + n(u')}f{E)-n(u'){l+n{u)}f(E')] , (1) 

where a is the fine structure constant, k = (uj,k) and 
k' — (w', k') are the initial and final CMB photon momenta, 
respectively, and p = (E,p) and p = (E',p') are the mo- 
menta for electrons. In Eq. CD, n{tj) and f(E) denote the 
distribution functions for the CMB photons and electrons in 
the CG, respectively, and X is the invariant transition prob- 
ability of the Compton scattering, which will be discussed 
in Section 3. In the present paper, we use the natural unit 
h — c = 1, unless otherwise stated explicitly. 

Assuming that the initial CMB photon distri- 
bution is spatially homogeneous, the space deriva- 
tive term in t he left hand side of E q . CD can 
be dropped. In Itoh. Kohyama and Nozawal (|l998T l and 
iNozawa. Itoh and Kohyama ( 19981 ). the Boltzmann equa- 



tion of Eq. CD was calculated in the CMB frame under this 
assumption. 

On the other hand, the left hand side of Eq. CD can be 
rewritten with variables in the CG frame as follows: 



dn(uj) ? _ , . dn c (uj c 



dt 



dt c 



+ k c ■ V c n c (id c 



(2) 



where k c = (tj c , k c ) denotes the momenta of the initial CMB 
photons in the CG frame, and n c (u) c ) is the CMB photon 
distribution function in the CG frame. Assuming that the 
initial CMB photon distributions are spatially homogeneous 
in both frames, one obtains 



dn(u>) 
' dt 



dn c (uj c ) 
dt c 



(3) 



Let us consider that the CG is moving with a velocity j3 c 
with respect to the CMB frame. Then, the photon energies 
are related with each other by the Lorentz transformation: 



^c = 7c (1 - Pcllc) LU , 



where 7 C = — /3i 

[x c = j3 C 'k 



and 



(4) 



(5) 



is the cosine between the unit vectors j3 c and k defined in 
the CMB frame. Inserting Eq. @ into Eq. C3), one has 



dn(x) . dn c {x c ) 

- 7c (1 - Pc/J-c) 



dt 



dt c 



(6) 



where we introduced new variables x c = Wc/fcalcMB and 
x = Lu/ksTcMB- Equation © gives the relation for the 
rate equations between the CMB and CG fra mes. Note that 
Eq. C51) agrees with the formula described in IChluba et al.l 

On the other hand, one can reverse the relation of 
Eq. by using the inverse Lorentz transformation: 



W = 7c (1 + /3cMc)^c , 

where 

t^c = Pc ' kc 



(7) 



(8) 



is the cosine defined in the CG frame. Instead of Eq. (|6]), 
one obtains 



dn c (x c ) , c , dn{x) 

= 7c (1 + /?cMc) 



d.t c 



dt 



(9) 



As far as the initial CMB photon distributions are spatially 
homogenous in both frames, Eq. (O and Eq. Q are equiv- 
alent. 

We now rewrite Eq. (|6]) with new variables r and r c as 
follows: 



dr c 



dn(x) dn c (x c ) 
- (1 PcM 

r = N e OT t , 

T c = NcCTT tc , 



(10) 

(11) 
(12) 

where ot = 8na 2 /(3m 2 ) is the Thomson scattering cross 
section, and N e and Nj; are the electron number densities 
in the CMB and CG frames, respectively, which are related 
with each other by 



N e = 7ciVe 



(13) 
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In Appendix A, we study relations of the electron number 
densities among different Lorentz frames, where we give the 
explicit derivation of Eq. (JTSJ). 

Finally, we extend the relation of Eq. ([2]) to the OBS 
frame moving with a velocity fi with respect to the CMB 
frame. Under the assumption of the spatial homogeneity for 
the initial CMB photon distributions, one has 



dn (x ) dn(x) 
— ^ = 7o (1 + PofJ-o) 



dt 



(14) 
(15) 



where j a — l/yl — Po, x a = cjo/fcsTcMB, and (uj ,k ) de- 
notes the four momenta of the initial CMB photons in the 
OBS frame. Inserting Eq. © into Eq. (fT4")l . one has 



dn (x ) . q. , dn c (x c 

~ Jojc (1 + PoUo) (1 - Pcflc) 



dt„ 



dtc 



(16) 



In Eq. (|16|) , k of fi c = Pc ■ k has to be transformed to the 
variables in the OBS frame. IChallinor and Leeuwer] (|2002r i 
discussed the transformation law for the photon propagation 
directions betwe en two different Lorent z fram es. One can 
reverse Eq. (2) of lChallinor and Leeuwer] (2002) and obtains 
as follows: 



fi ■ v — 3 \ . n' — n' ■ vv 
1 v + 



7 (1 — h' ■ v) ' 



(17) 



Inserting n = — k, n = — k and v — B into Eq. (I17|l . one 
obtains 



Po • k -\- do i S _i_ ° @° ' ^° @° 

Po T 



k - 

l+0 o -ko 
Thus, one has 

(l + A^Hi-A^c) 



7o (l + (3 B ■ k ) ' 



(18) 



= (1 + /3a ■ ko) 



1 



Po ■ k a + Bo 



Be- Po 



l + Po-k 

fic ' k Q $o ' k Q $c ' fto 



7o (l +P a ■ k ) 

« (l - Pa ■ Pc) [1 - {Pc£ ~ PoH°a)] , 
where 

fj, c — f3 c ' k Q • 



(19) 



(20) 



In deriving the last equality in Eq. (|19[1 . 0(p c Po, PoPc) terms 
were neglected. One has 

—^ t = Jo 7c (1 - Pa ■ Po) (1 - PcfJ-c + PoUo) -(21) 



Finally, we rewrite Eq. ()21[) with a new variable r as 
follows: 



drio(xo) „ ^ dn c (x c ) 

j — l- 1 - Pc^c + Pof^o) 

(LTo 

t = N°a T t , 



dr c 



(22) 
(23) 



where 7V° is the electron number density in the OSB frame, 
which is related to N£ by 

K = 7c 7c (1 - Pa ■ Pc) K . (24) 

In Appendix A, we give the explicit derivation of Eq. (|24[) , 



Thus, Eqs. ()10|) and ()22[) summarize the formal relations 
for the rate equations among the CMB, CG and OBS frames. 
In Section 3, we study the rate equation in the CG frame. 
Then we transform the rate equation to the CMB frame with 
Eq. (fTOl) in Section 4 and to the OBS frame with Eq. 
in Section 5. 



3 CALCULATIONS IN THE CG FRAME 
3.1 Rate equation 

Let us first define the kinematics in the CG frame. The four 
momenta of the initial and final electrons are p c = (E c ,p c ) 
and p' c = (E'cPc), respectively, and P = p c /E c is the veloc- 
ity of the initial electron. The four momenta of the initial 
and final photons are k c = (oj c ,k c ) and k' c = (cj' c ,k' ), re- 
spectively. We introduce the direction cosines of the angle 
between the initial electron and photons as follows: 

^ = p ■ k c = cos 9 kc , (25) 

fi = p ■ k c — cos 0^ . (26) 

According to iNozawa and Kohvamal l|2009h . the rate equa- 
tion in the CG frame can be written in the Thomson ap- 
proximation as follows: 



dn c (x c 
dr c 



64tt 2 



dpcP c Pe,c(E c ) j dQ Pc j dQ k / c -^ 



- ) X A [n c (x' ) - n c (x c )] , (27) 



1 — P(JL \ Ul, 



<£c_ _ 1 - Pii 



Pn> ' 



(28) 



where 7 = 1/ yjl — P 2 , and the electron distribution func- 
tion is normalized by dp c PcPe,c{E c ) = 1. In Eq. (|27[) . 
Xa is given by 



X A 



where 



2(1 - cose e ) (i - cose c ) 2 

7 2 (i - Pv)(i - M) + 7 4 (1 - M 2 (i - Pv'Y ' { ' 



s9 c = + y/l - fl 2 y/l - fl' 2 cos(0 fcc - (j> k i) 



(30) 



is the cosine between k c and k' c . It should be noted that 
the Thomson approximation used in deriving Eqs. (I27|l and 
l|28p is extremely good approximation for the CMB photon 
energies. 

Now, we proceed the calculation by expanding the pho- 
ton distribution functions n c (x c ) and n c (x' c ) in power series 
of p c . They are Lorentz invariant, namely, 



n c (x c ) = n(x) , 
n c (x' c ) — n(x') , 
where 

X = 7c(l + PcfJ.c)Xc , 
X = 7c(l + PcH'c)Xc , 

[il = p c ■ k c = cos 8c , 
He' = Pc • k'c — cos 9c' . 



(31) 
(32) 

(33) 
(34) 
(35) 
(36) 
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The angles fj, c c and /i c are related with each other by the 
Lorentz transformation: 

= -, o ' ( 37 ) 

1 — (j c fl c 

Inserting Eqs. J33J and ((34j into Eqs. J31]) and (|32)l . respec- 
tively, and expanding up to O{0£) terms under the assump- 
tion /3 C <C 1, one obtains as follows: 

n c (x c ) = n(x c ) + ^/3 2 D Xc (D Xc + 2) n(x c ) 



+ \plP2{£) D Xc {D Xc -1) n(a; c 



(38) 



n c (sc) = n(x' c ) + -PcD x , c (D x > e + 2) n(x' c ) 
+P c Pi(jie) D x ,n{x' c ) 

+~ P >lP 2 { l £)D < {D x , c - l)n{x' c ) , (39) 

where Pi{^l c c ) and Pe(fic') are the Legendre polynomials of 
the £-th order, and 

d 



D 



dz 



(40) 



is the Lorentz invariant operator. 

Choosing the direction of f3 c along 2-axis, one can reex- 
press Pi(fj,c') as follows: 

m 

= P«(cos9 c )PtQit), (41) 

where 3> c and $i are the azimuthal angles of fc c and fc c , 
respectively. In deriving Eq. (|41|) . only m — term survived, 
because the photon distribution of k c is axially symmetric 
for the z-axis. 

Now, we introduce a new variable s c by 

x' c = e Sa x c . (42) 

Inserting Eqs. l[3"gj> . (O and (gj) into Eq. l(2"7)). one finally 
obtains the rate equation in the CG frame as follows: 



dn c (x c 
dr c 



ds c P , c (s c ) [n(e Bc x c ) - n(x c )] 



+ I ds c P , c {s c )D Xc (D Xc +2)[n(e 3c x c )-n{a 



+ /?cA(Mc 



where 



dsc Pi,c(s c ) D Xc n(e 3a x c ) - D Xc n(x c 



ds c P 2 ,c(s c ) D Xc (D Xc - 1) n(e Sc x c 



-D Xo (D Xc - l)n(x c ) 



(43) 



ds c Pf,c(s c )n(e Sc 3; c ) 



64tt 2 



dp c p 2 p e , c (E c ) J dQ, Vc J dUy^ 

x /a 2 X* Pi{cosQ c ) n{x' c ) . 
(1 - Pfi'y 



(44) 



In deriving Eq. (|43|) , we used the normalization for Po,c(s c ), 
rfs c Po : c(s c ) = 1. (45) 

In Appendix B, we study on the normalizations for Pe, c (s c ). 

As we will show explicitly in Section 6, the redistribu- 
tion function Pe, c (s c ) satisfies the following useful relation: 



Pe,c(s c ) = e 3s " Pi,c(s c ) 



(46) 



for £—0, 1 and 2. Because of the relation, one can use the 
same rate equation of Eq. (|43l) for the spectral intensity 
function I(x c ) = x 3 n c (x c ) /2tt 2 except for replacing e s " by 
e~ Sc and D Xc by D Xc — 3. Namely, one has as follows: 



dl c (a 



dr c 
6 



+ /3 c Pi(Mc 



+ -P 2 C P2^ C C 



ds c P , c (s c ) [j(e s "x c ) - I(x c )\ 

ds c P , c {s c ) (D Xc -3){D Xc -1) 

x [I(e-'x c ) - I(x a )] 
ds c Pi,o(s ) {D Xe - 3) I(e~ a "x c ) 

-(D Xc -3)I(x c 

-OO 

da c P 2 , c (sc) (D Xc -3)(D Xc -4) 

SO 

xl(e- s "x c ) - (D Xc - 3)(D Xc - 4) I(x c ) 



(47) 



It should be emphasized that the property of Eq. (|46[) for 
£=0, 1 and 2 are very remarkable features of the CG frame. 



3.2 Photon number conservation 

In the present subsection, we verify the photon number con- 
servation of Eq. ((43} by integrating over the phase space 
volume element d s x c for the initial photon. One has as fol- 
lows: 



l3 dn c (x c 
d x c 



6 



dr c 

ds c Po,c{sc) e~ 3s " J d 3 x' c n(x' c ) — J d 3 x c n(x c ) 

/ + oo n 
ds c Po, c (s c ) e" 3sc / d 3 x' c D x , c (D x , c +2)n(x' c ) 
- OO J 

d 3 x c D Xc (D Xc +2) 
+ .i.2n I dx c x 2 J dfi c c Pi(a*c) 

ds c Pl jC (s c ) D Xc n(e Sc x c ) - D Xc n(x c ) 



+ 3^ 2lV J dXcX c J d if c P2{n c c 



ds c P 2 , c {s c ) D Xa (D Xc - l)n(e 8c x c ) 
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-D Xa (D Xc - l)n(xc) 



(48) 



where we introduced x' c = e 3a x c . Note also that the z- 
direction was chosen along /3 C direction for the integration 
over the solid angles. Using Eqs. ((45} and (|46|l , the first two 
terms in Eq. (|48p cancel out with each other. The third and 
fourth terms also cancel out with each other by the same 
reason. The other terms become zeros by the ^ integrals. 
Thus, the photon number is conserved in the CG frame, 
namely 



, 3 dn c (x c 
a x c 



dr c 



= 0. 



(49) 



3.3 Operator representation and the 

Fokker-Planck expansion approximation 

In lNozawa. Kohvama and Itohl (|2009l ) . a formal solution was 
shown for the rate equation of the leading-order SZ effect by 
introducing an operator 0(D)= f_°° dsPi(s)e sD — 1. In the 
present subsection, we extend the operator representation 
formalism to the rate equation of Eq. (|43[> . Let us first in- 
troduce an operator Oe(D z ) by 



Ot(Dz) 



ds c Pi, c (s c ) e" 



(50) 



Then, Eq. ((43} can be rewritten as follows 

dn c (x c ) 



dr c 



O (D 
+ /3 c P 1 (p' : c )0 1 (D 



+ i/? c 2 O {D Xc ) D x 



D Xc + 2) 



+ -P 2 c P2(^)0 2 (D Xc )D Xc (D x 



n(x c ) . (51) 



We now derive the expressions for the opera- 
tors Ot(D z ) in the Fokker-Planck expansion approxima- 
tion. The expression for Oq(D z ) was already shown in 
iNozawa. Kohvama and Itohl (2009). We recall the results 
here so that the present paper is self-contained. They are 
as follows: 



O (D M 


) = Oo,n9" , 
n = 


(52) 


00,0 = 


o, 


(53) 


00,1 = 




(54) 


00,2 = 


3 . 7 a2 

-Io A * + io A - 


(55) 


00,3 = 


120 15 z 30 z ' 


(56) 


00,4 = 


23 . 23 . 2 431 * 3 16 A 4 

— A z H At At H At , 

56 21 z 420 105 ' 


(57) 


00,5 = 


1457 . 309 . 2 839 A 3 74 . 4 

A z At H A, A, + 

896 140 336 105 


^A«,(58) 



D z (D z +3) 



J_d_ 

z 2 dz 



(59) 



where 8 e = ^bIcmb/™ and 0(6®) terms were neglected. It 
should be emphasized that the photon number conservation 
is automatically guaranteed for the operator Oq(D z ) in the 



Fokker-Planck expansion approximation, because it can be 
written as a function of A z alone. 

Similarly, one obtains the expressions for Oi(D z ) and 
2 {D Z ) up to 0{6l) term as follows: 



(i) for i = 1 

5 

0l(Dz) = ^01, 



01.0 = -1, 

01.1 = -|-|a z , 

5 5 

12 2 2 

01.2 = -- + -A, - -A\ , 

407 233 . 34 A 2 17 A 3 
= 140 _ 140 + 35 _ 70 



28 84 



21 



420 



210 



_ 23363 160165. , 3547 . 2 2143 . 3 

^1,5 = —T77, + r^rA, 



448 4032 



252 



560 



+— A* - —Al 
45 126 



(ii) for £ = 2 



2 {D Z ) =^02,nC, 





71= 







02,0 = 


9 

10 ' 






02,1 = 


3 

-5 + 






02,2 = 


183 
70 


|a, + 1a? 




02,3 = 


429 

_ lo~ 


1535 . 

H A z - 

336 


14 A 2 23 . 3 
— H At , 

15 210 ' 


02,4 = 


2559 
56 


8287 . 

A z + 

336 


17 A 2 341 . 3 

— At At + 

3 420 



35 



26961 704183 . 70843 . 2 43901 . 3 

^2,5 = 7^7 1 ^ n „„ A z n ^ ^ A z + _ 



128 5376 



2100 



8400 



18 A 4 4 a5 

- — At H A z 

35 175 



(60) 

(61) 
(62) 

(63) 
(64) 



363 703 . 67 A 2 409 A 3 23 . 4 
1A = -— + —A z ~-At + —Ai-—A 4 z , (65) 



(66) 

(67) 

(68) 
(69) 
(70) 

(71) 

(72) 

(73) 



Note that unlike the 00(^2) operator, the Oi(D z ) and 
02 (-Dz) operators contain nonzero constant terms. There- 
fore, the photon number conservations are not automatically 
guaranteed for these operators. 



4 TRANSFORMATION TO THE CMB FRAME 
4.1 Rate equations 

In the present section, we rewrite Eq. (|43jl in terms of vari- 
ables in the CMB frame. Inserting the result into Eq. (|10[) , 
one obtains the expression for the rate equation in the CMB 
frame. 

First, we expand x c up to O(0 2 ) terms under the as- 
sumption /? c < 1 as follows: 
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(74) 



1 , 

X c = X - P c fl c X + ~P C X , 



Then, one has 

n(x c ) = n(x) + ^p 2 D x (D x + 2) n(x) 
+PcPi(p c ) D x n(x) 



n{e s "x c ) = n{e s <=x) + ~p 2 D x (D x + 2) n{e s "x) 

+/3 c P 1 {fi c )D x n( e s,: x) 
1 



(75) 



+-PiP 2 ^c)D x (D x - l)n{e"x). (76) 



Equation (|37|l can be also expanded as 

2 
3 



A = + % {ft(/io) ~ 1} Pc + 0(Pl 



(77) 



Inserting Eqs. I|75 |l -(|77 [l into Eq. (|43|l. one has 

dn c (x,^ c ) 



3 

-PcPx(»c 
+ \plP2{»c) 



ds P m (s) [n(e s x) — n(x)] 
dsP d (s)D x (D x + 2) n(e s x) 

+ oo 



dsP d (s) Arn(e s a;) 

-oo 

dsP q (s)D x (D x - l)n(e s x),(78) 



where we introduced 
P m (s) = Po,c(s), 
Pd(s) = Po,c(s) - Pl,c(s) , 

=P2,c(s) -2Pi,c(s) + P ,c(s). 



(79) 
(80) 
(81) 

In deriving Eq. (|78[) . we replaced the integral variable s c by 
a new variable s. 

Before proceed to the final step in deriving the rate 
equation in the CMB frame, it should be noted as follows. 
As shown in Section 3.2, Eq. (|43|l satisfied the photon num- 
ber conser vation. However, Eq . (|78p which corresponds to 
Eq. (25) of lChluba et~all l|2012t ) contains the terms violating 
the photon number conservation in the p 2 term. This can be 
checked immediately with the explicit forms for the Oi(D z ) 
operator of Eqs. (|60|) - (|66l) . This illustrates that the pho- 
ton number conservation is representation dependent. On 
the other hand, we will show in the next subsection that 
Eq. (|10f) conserves the photon number. 

Inserting Eq. (|78|) into Eq. (|10[) . one finally obtains the 
rate equation in the CMB frame as follows: 



dn(x) 



dr 



+ ~Pc 

3 



ds P m (s) [n(e 3 x) — n(x)] 
ds P d (s) D X (D X + 3) n(e s x) 



PoPl{Hc) 



+ 



dsPd{s)D x n(e a x) 

O 

ds Pm(s) {n(e s x) — n(x)} 



+ \PIP2{^) 



+ 2 



dsP q (s)D x (D x - l)n{e 



ds Pd(s) D x n(e s x) 



(82) 



In the next subsection, we compare the present result with 
previous works. As shown in Eq. ((79}, the expression of the 
leading-order thermal SZ effect is common for both frames. 

Finally, it should be noted that one can use Eq. (|82[) for 
the spectral intensity function I(x) = x 3 n(x)/2ir 2 except 
for replacing e s by e~ s and D x by D x — 3. This can be 
immediately verified with Eq. (I46|) . 



4.2 Photon number conservation 

In the present subsection, the photon number conservation 
of Eq. (|82p is checked by integrating over the phase space 
volume element d 3 x for the initial photon. One has as fol- 
lows: 

d\^fil 
dr 

ds P m (s)e _3s J d 3 x' n(x') — J d 3 xn(x) 
+ \pl [ dsP d {s)e- 3a ( d 3 x' D x i (D x i + 3) n(x') 



ds Pd(s) D x n(e 3 x) 



— P c 2ty J dxx J dfi c Pi(fj, c ) 

— / dsP m (s) {n(e"x) — n(x)} 



1 2 / 2 / 

+ ~P C 2% J dxx J dfi c P2(nc) 



dsP q (s) D X (D X ~ 1) n(e B x) 



+ 2 dsP d {s)D x n{ex) , (83) 



where we introduced x' = e s x, and therefore, d 3 x = 
e~ 3s d 3 x'. Note also that the ^-direction was chosen along 
Pc direction for the integration over the solid angles. 

Using Eqs. (gSJ and (gBJ , the first two terms in Eq. (fg3")l 
cancel out with each other. The third term becomes zero, 
because 



/ dx x' 2 D x i (D x i + 3) n(x') 
Jo 



fo ^ dx ' 


,4 dn(x') 


, 4 dn(x') 


OO 




X dx ' 


X dx' 






= 0. 



(84) 



Other terms become zeros by the [i c integrals. Thus, the 
photon number is conserved in the CMB frame, namely 



dr 



4.3 Comparison with previous works 



(85) 



In orde r to compare the present work wit h previous works 
such as lNozawa. Itoh and Kohvamal (j 1998h (denoted as NIK 



© 0000 RAS, MNRAS 000, [T^?? 



Analytical studies on the SZ effect 7 



as under) and IChluba et all fcoij ), we first reexpress the 
rate equation of Eq. (|82|) in the operator representation with 
Eq. (|50[) . One has as follows: 



dn(x) 
dr 



O (D x ) 

+ ^/3 2 {O (D x ) - Oi(D x )} D x (D x +3) 

- /3 c Pi(fi, c ) {O {D x )(D x + 1) - Oi(D x )D x } 

+ \plP2{Hc) {0 2 (D x )D x (D x - l)-20i(D x )Dl 



+O (D x )D x (D x + 1)} 



n(x) . (86) 



On the other hand, one can also express the rate equa- 
tion of NIK (1998) in the operator representation. One has 
as follows: 



dn(x) 
dr 



1 

7c 
1 



dn(x) 



dr c 



Oo{D x ) + i/3 c 2 Onik,3(C 
- PcPiific) OmK,i(D x ) 

+ -^P 2 (/i c )0NIK,2(A c 



n(x) . 



(87) 



In Eq. ()87|) . the factor 1/7 C was multiplied to the relate 
equation of NKI (1998) in order to be consistent with the 
definition of the rate equation of Eq. (1861) . The reason for 
the factor is because N£ was used in NIK (1998) instead 
of N e for the electron number density in the CMB frame. 
The factor is consi stent with the comment indicated by 
IChluba et al.l {2012). The operators OmK,e(D z ) are calcu- 
lated by Eqs. (14)-(17) of NIK (1998). They are 



Onik,i(-D z ) = ^ ' Oi, n 6" ■ 
D z 



01,0 

01.1 = 

01.2 = 

01.3 = 

01, 



5 5 5 

7 _, 39 ^ 2 67 ^3 73 ^ 4 

D z H DZ H Di H D z 

10 10 5 10 



11 

Io J 



D z 



1031 
280 
3749 
210 

795 



373^ 2 

84 z 



+ 



56 
2979 



381^2 

28 z 



840 
_64 

210 105 
7961 



9607^ + 2939^ 



210 



(89) 
(90) 

(91) 
(92) 



840 



1469 6 489 7 
84 2 35 



i 



140 

0NIK,2(^ Z ) = \ 02, n^? . 

°><° = To D » + To D - 

• 5 10 5 10 



17873 ^4 
420 z 



(94) 

(95) 
(96) 



0NIK,3(-Dz) 



E 



(97) 



03,0=^(^+3), 

/ 1 Q 91 7 

03,i = D Z (D Z + 3) (- + -D z + -D\ 



(98) 
(99) 



Comparing Eqs. (JHU+flMJ with Eqs. fl52])-(E3J , one 



finds 

0nik,i(D,) = O (D z ){D z + 1) - O x {D z )D z 



(100) 



0NiK, 2 (£g = 2 {D z )D z (D z - 1) - 2 Oi (£>,)££ 

+ 00(^)^(^ + 1). (101) 

Equations (|100|l and (|101[1 show that two calculations agree 
for the p c term and /3 2 P2(Mc) term. On the other hand, as 
for the 0nik,3(-Dz) operator, there is a difference in Eq. ()99|) 
by amount of 3/2D z (D z + 3). However, this term exactly 
cancels with the /3 2 term in the l/7 c 0o(Az) term, namely, 

- i/3 c 2 0o(D z ) = ^D z (D z +3)e e + 0(6 2 e ) . (102) 

Therefore, the two calculations agree completely. This clar- 
ifies the relation of the calculations between the CMB and 
CG frames. Namely, as far as the CMB photon distribu- 
tions are spatially homogeneous in both of the CMB and 
CG frames, the two calculations are equivalent. 



5 TRANSFORMATION TO THE GENERAL 
OBSERVER'S FRAME 

Let us now consider that the OBS frame is moving with a 
velocity f3 with respect to the CMB frame. The variables in 
the OBS frame are defined in Section 2. The photon energies 
x and x are related by the Lorentz transformation: 



X = Jo (1 + /3ofJ.°)Xo ; 



(103) 



where [i° is given by Eq. Then, one expands n(x) and 

n(e s x) up to 0(/9 2 ) terms under the assumption /3 -C 1: 

n(x) = n(x ) + ^P 2 D XO {D Xo + 2) n{x„) 
+/9 Pi(/i°) D Xo n(x ) 

+ i/3 2 f 2 (Mo) D Xo (D Xo -l)n( Xo ), (104) 



n(e"x) = n(e a x ) + -/3 D Xo (D Xo + 2) n(e s x ) 
6 

+f3 P 1 {^ )D Xo n{e a x ) 

+ \plP2{n° ) D Xo (D Xo - l)n{e 8 x ) . (105) 

Similarly, /3 c Pi(fMc) and /3 2 P2(/i c ) can be also transformed 
with Eq. (|17(l as follows: 

PcPl(Vc) = PcPlQg) +Po-Pc-0o- kjc ■ k , (106) 

/3 2 P 2 (Mc) = PlP2{^) , (107) 



where is given by Eq. (|20l) . 

We finally obtain the expression for Eq. (|22|l with 
Eqs. (|T04|) — (fTUT Tl as follows: 



drio(xo) 
dr 



dn {x ) 
dr 



+ 



dn D (a;o) 
dr 



+ 



drioixo) 
dr 



,(108) 
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where 



dn {x ) 



dr a 



dsP m (s) \n(e s x ) — n(x )] 



+ -/?./ ds P d (s) D Xo (D Xo +3)n(e s x ) 



+ ^ 2 C Pl(H°c) 



dsP d {s) D Xo n(e s x a ) 
+ / dsP m (s) {n(e s x ) —n(x )} 

SO 

ds P q (s) D Xo (D Xo - l)n{e"x 



+ 2 dsP d {s)D Xo n(e s x ) 



(109) 



dn {x ) 



dr a 



= -Pl / ds P m {s) D Xo (D Xo +4) [n(e a x ) - n{x )] 

J — oo 

+ /3 o Pi0°) / dsP m (s) (D Xo + l) [n(e s x )-n(x )] 



+ ^oP2(fi°„) I dsP m (s)D Xo {D Xo + l) 



x [n(e s x ) - n(x )] , (110) 



dn {x ) 



dr 
1 



A, /3c 



+ 



dsP d (s) D Xo (D Xo + 3) n{e a x ) 



dsP m (s)D Xo {n(e"x ) — n{x )} 



Po ' hope ' ko ' ^ c 



ds P d (s) D Xo n(e s x ) 



+ 



dsP m {s)D Xo {n{e a x ) —n(x )} 



(111) 



The photon number conservation is manifest for the 
rate equation of Eq. ()108p by the following reason. Equa- 
tion (I109[) which is identical to Eq. (|82[) satisfies the photon 
number conservation. The photon number conservation of 
Eq. (fTT0|) is shown with Eqs. (g5j) and ([46]). As for Eq. (fTTTj) . 
the first two terms in the [] brackets are zeros by the same 
reasons discussed in Section 4. The last two terms become 
zeros by the integral over the solid angles, because the terms 
in the ( ) brackets are expressed by 



Po ' koPc ' k Q ' @ c 



1 /8tt 
5 V 15 



+ 2 



J2 (-l) m [/3o x Y 2 ,- m (k ) 



(112) 



Therefore, the photon number is conserved for Eq. (|108[) . 
namely 



d s Xo dn (x ) =Q 
dr 



(113) 



Let us now study the expression of Eq. (|108[) in the limit 
Po — Pc- It is nontrivial whether the obtained expression has 
the correct form in the limit p o = p c . In this limit, one can 
show that fi° — fi° — n c c and N° = A^^ (shown in Appendix 
A). Then, it is straightforward to show that 



dn (x ) dn c (x c ) 



dr 



dr c 



(114) 



Therefore, Eq. (|108[) is correctly reduced to Eq. (|43|) in the 
limit of Po —¥ p c as it should. 

Before cl osing the prese n t sect ion, it should be remarked 
as follows. In IChluba et all (|2012T ). it was suggested to use 
/x c instead of p° in describing the rate equation of Eq. (|108[) . 
These angles are related with Eq. (|18p by 



fl° = He + Po A*c Mo - Po ■ Pc + O(Po) ■ 



(115) 



Inserting Eq. (|115|l into Eq. (|108p . one finds the following 
modifications. As for Eq. (|109[) . Pe(fJ-c) should be simply 
replaced by Pi(fj, c ) for £ = 0, 1, 2. There are no modifications 
to Eq. <)11U| I. On the other hand, Eq. (|llip is modified as 
follows: 



dn (x ) 



dr 
1 



PoP, 
Po ■ Pc 

+ 



dsP d (s)D Xa (D Xa +l)n(e 8 Xo) 



dsP m (s) (D Xo - 2) {n(e s x ) - n(x )} 
[Po -kPc-k- ~p o ■ Pc^j 

ds P d {s) D Xo (D Xo + l)n(e 3 x ) 

+ 1 dsP m {s){D Xo + 1) {n(e s x ) - n{x )} .(116) 



It is clear that the first term in the Po ■ Pc terms does not 
satisfy the photon number conservation unless p o ■ Pc = 0, 
which again illustrates that the photon number conservation 
is representation dependent. 



6 CALCULATION OF Pt,c(So,P) 

6.1 Analytic expression for P^, c (s c , /?) 

For the practical calculation of the rate equation, Eq. Q43p 
can be solved by calculating 5-dimensional integrals of 
Eq. (|44p . For the present computer technologies, the direct 
numerical integration can be done within a reasonable CPU 
time. However, it is extremely important to derive analytic 
expressions not only to reduce the CPU time of the com- 
puter, but also to investigate various properties of functions 
in the equation. The obtained analytic expression will be 
extremely useful to analyze observational data. 

The purpose of the present subsection is to calculate the 
redistribution function Pe t c(s c ) in Eq. (|44p . The calculation 
is simplified in the electron rest (ER) frame. According to 
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iNozawa and Kohvamal ((2009), this can be done by the fol- 
lowing Lorentz transformations for the photon angles from 
the CG frame to the ER frame: 

-fio + P 



" 1 - /W 

/ _ -Mo + P 
11 1-/>A*&' 
where 

= $ ■ k , 

Mo = P ■ k' ■ 



(117) 
(118) 

(119) 
(120) 

In Eqs. (|119[) and ()120[) . ko and k' are unit vectors of the 
initial and final photons in the ER frame, respectively. Then, 
Eq. (|29|) is simplified as 

Xa,o = | Po(cose ) + | P 2 (cose ) , (121) 
where 

cosOo = jUo/u'o + \A - A*o \/ 1 - Mo 2 cos (0fcc - ^kj.) (122) 

is the cosine between fco and fco- In Eq. (|121[) . the first term 
corresponds to the isotropic scattering and the second term 
gives the anisotropic scattering. 

Similarly, Pe(cosO c ) can be also reexpressed with 
Pe(cosOo) as follows: 



Po(cos6 c ) = P (cose ) = 1, 
Pi(cosS c ) = 1 — a + aPi(cosOo) , 

P 2 (cos9 c ) = 1 - 3a + 2a 2 + 3a(l - a) Pi(cos9 ) 



+0 P2(cOS0o) . 



where 



7»(l-/9/io)(l-/9/4) ■ 
Inserting Eqs. (frT?) and (flTS)) into Eq. pS). one has 

e "c _ ^4 _ 1 - PfJ-'o 

LJc 1 - 0fM> ' 

which implies 
Mo = ^[l-e ac (l-/?Mo)], 

rf^o = -4 (1 - Pt*o) e Sc ds c ■ 

Thus, fi' is expressed with fj,o and s c . 

One can then reexpress Pe,c(s c ) in Eq. (|44[) as follows: 

Pf, c (s c )= / d/3/3 2 7 5 p e , c (/3)P f , c (s c ,/3), (130) 

1 - e~ |3cl 
1 + e l Sc l 

where pe,c(P) = m 3 p e ,c(E c ), m is the electron rest mass, and 



(123) 
(124) 

(125) 
(126) 

(127) 

(128) 
(129) 



3e s 



64tt 2 /37 4 



C2(«c) 



d/io 



(l-/?Mo) 2 



'Ml( s c) 

f27T / i 27T 

x / d<£ fcc / d<f> k , X A ,o P<(cos9 c ) , (132) 



/o ^0 
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/ \ i 1 for s c < , . 

W*>-\ [l-e- s <=(l + /3)]//3 fors c >0 ' {M) 



= < f-^-w* ; :;: <o , (I34) 

Before proceed to the practical calculation, it should 
be emphasized as follows. In Eq. (|132[1 . photon zenith an- 
gles (fio and fi' ) are described in the ER frame. On the other 
hand, energies (cj c , u)' c and P) and azimuthal angles (<f>k c and 
4>y ) are left in the CG frame. Therefore, the present coor- 
dinate frame is a hybrid system of the CG and ER frames. 
When the triple integrals are completed, however, the ob- 
tained Pe,c(s c ,P) becomes the function of only variables in 
the CG frame. 

Inserting Eq. (|12ip into Eq. (|132p . one has 



PeAsc P) = PeoAsc P) + P«2,c(sc, P) ■ 



(135) 



PtoA s dP) = 



3e s 



i6/3 y 



M2(s c ) 



dfi 



' t ho , (136) 



Cl(sc) 



(1 - /? M o) 2 3 



3e s 



PaA*°>0) 16/ 3 7 4 
where 



Iff — 



M2(sc) 



(2tt)» 



^° (i-L) 2 i J "' (137) 



P*(cos9 c ) P f -(cos6o) .(138) 



Equation (|138[l can be integrated analytically, one obtains 
as follows: 



(i) for i = 
loo = 1 , 

Jo2 = P2(a*o)P2(Mo) , 

(ii) for £ = 1 

/io = 1 + a (-1 + MoMo) 



(139) 
(140) 

(141) 



f 2 

li2 = P2(mo)P2(mo) + « I g MoMo ~ P2(Mo)Pi(Ato) 

+ ^P 3 (Mo)P 3 (Mo)} ,(142) 

(hi) for £ = 2 

/20 = -2 + 3/io + a 2 {2-3^oMo + P2(Mo)P2(Ato)} , (143) 



/22 = -2P 2 ( MO )P 2 (Mo) +3/i2 

1 6 , 16 



+a 2 {- - - MoAto + — P 2 (mo)P2(mo) 

-| P 3 (mo)P 3 (mo) + § P 4 (/io)P4(Mo)} , (144) 



where a and /i are given by Eqs. (|126|) and (|128|) . respec- 
tively. 

Our final task in the present subsection is to derive 
analytic expression for P U i , c (s c , /3). After lengthy but 
straightforward calculation, one finally obtains the follow- 
ing results: 
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Poo,c{s c ,/3) = 



28 2 j 2 



8 (cosh— — — sinh-^J- 

' \ 2 8 2 



(145) 



P02, C (S C ,£) = -Poo, c (Sc,/3) 



"4/3 2 7 2 



3 / 2 S C 1 . 2 |«c 

—— - cosh smh J — 

8-y 2 \ 2 B 2 2 

| 3cosh-^ — — sinh^-! 



9 / s c 1 |s c | \ 2 s c 

H — = — 7 cosh smh- — 1 cosh — 

/3 3 7 4 V 2/3 2 1 2 



' (3-/3 2 )(A^-| Sc |)cosh| 



where 
A/3 = In 

cosh a = 
sinhz = 



2/3 4 7 

1 + /? ' 

e z + e 
2~ 

e z -e 



(146) 

(147) 

(148) 
(149) 



Note that Eq. (|145[) is the familiar expression in the isotropic 
scattering approximation. On the other hand, Eq. (|146[) 
gives the anisotropic scattering contribution. 
Similarly, one has as follows: 

(ii) for I = 1 

Ao,c(s c ,/3) 

= 6 „ —8 ( cosh — — — sinh^- 
2^ 2 7 2 3 I 2 8 2 



If ,S C 1 \S C \ \ . 2 \Sc\ 

_f CO sh T --smh— jsmh — 



(150) 



Pl2,c(Sc,8) = -P W , c (Sc,B) 



a 3s c /2 



3 / 2 Sc 1 . ,2 |Se 

—— - cosh smh J — 

8-y 2 \ 2 B 2 2 

x | 5 cosh— — — sinh^l 
V 2 8 2 



{^i - ^ sinh ^ ) l 82cosh f 

25 ,8 C . ,\s c \ 20 . , 2 |s c | 

+_ cosh _ smh ___ smh _ 

75 / s c 1 . |s c |\ 2 s c 

H — = — =■ cosh smh- — - cosh — 

/3 5 7 6 ^2/3 2 J 2 

(A 8 - |s c |) (3 - /3 2 ) cosh^ 



2^ 4 7 2 
3 



i cosh- 



2/3 6 7 4 (^-M){3(5-/? 2 ) 

(5 -3/? 2 ) cosh 3 !} 



+2 i 



.(151) 



(iii) for £ = 2 

fWSc,/?) 

= 6 „ —8 | cosh— — i sinh-^- 

2^ 2 7 2 5 P I 2 ^ 2 



,s c 1 \s c 

h 2 v cosh y-^ sinh ^- 



cosh — smh 



H — -!— - [cosh— — i sinh-^J- ] sinh 4 -^-!- 

£ 3 7 4 \ 2 3 2 J 2 



2 

(152) 



P 2 2, c (s c ,/3) = -P 20 ,c(Sc,/3) 



3s c /2 



4/3 2 7 2 



07* 



cosh 



2 s c 



1 . , 2 |*c| 

^2 smh — 



x I 3 cosh— — — sinh-^- 

2 8 2 



I 3 fcosh 2 — - -Lsinh 2 li^O 

/? 3 7 4 V 2 P 2 / 

x f 94 cosh 3 £2 - ^ cosh 2 ^ sinhM 



2 

3 l s c 



54 ,! t .,j s c 18 . , 
--cosh- smh _ + _smh ^ 



, cosh— — — sinh— 

4/3 5 7 6 \ 2 8 2 



x ( 32 cosh 2 — — ^ cosh— sinh-^ 
2 8 2 2 



40 . ,2 MA . ,a|«c| 



15 



i 1 . , \Sc 

. cosh smh — 

4/3 5 7 6 \ 2 8 2 



x 211 cosh cosh — smh J — - 

1 2 /3 2 2 



+ 



2205 
4/3~V 
3 

" 2/3 4 7 2 
9 

" 2/3 6 7 4 



161 . ,2 |Se, , , ■ 

_ smh — jcosh 2 

, S c 1 . , |S C | \ 4 S C 

cosh sinh J — - cosh — 

2/3 2 ) 2 



(As - |s c |) (3 - /3 2 ) cosh^ 



(A/3 



! cosh- 



45 (X 



s c |)[3(5-/3 2 ) 

+2(5-3/3 2 )cosh 3 ^} 
5c |){3(7-/3 2 )cosh^ 

+4(7-3/3 2 )cosh 3 !} 



.(153) 



Thus, 5-dimensional integrals for the rate equations are re- 
duced to 2-dimensional integrals not only for the calculation 
in the CG frame but also in the CMB and OBS frames, which 
will reduce the CPU time of the numerical calculations dras- 
tically. 
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6.2 Properties of functions Pe,c(s c ,P) and Pe,c{s c ) 

In the present subsection, we study various properties of the 
redistribution functions Pe iC (s c ,fi) and Pe, c (s c ). As already 
indicated in Section 3, Pe,c(s c ,/3) and Pe, c (s c ) satisfy the 
following relations: 

Pt, B {a B , P) = e 33 " P e , c (-s c , 0) (154) 
and 

Pf, c (s c ) = e 3s <=P,, c (-s c ). (155) 

In Eqs. (fTio) , (fT46|) and (fT50)) - (fT53l , all functions in the 
] brackets are even functions of s c . Therefore, a multi- 
plicative function e 3s "^ 2 in P M / c (s c ,/3) gives the property 
of Eq. (|154l) . Equation (|155[) can be immediately obtained 
with Eq. 1(150) 1 . 

Now, we study the behavior of the function Pp jC (s c ,/3) 
at specific values of /3. First, one has at /} = /3 m i n , 

Pe,c( S c, Pmin) =0, (156) 

for £=0, 1 and 2. This can be easily verified with the follow- 
ing identity relations: 

cosh^ - J— sinhM = o (157) 

^ Pmin ^ 

and 

^-N = ' ( 158 ) 

Then, one has at /3=1 as follows: 

7 2 P , c ( Sc ,l) = 3 e 3, c /2 e -KI/2 ) (159) 

7 2 P 1 , c ( Sc ,l) = |e 3 -/ 2 ie- 3 l-l/ 2 , (160) 

7 2 ^, c ( Sc ,l) = ?e 3s = /2 i e - 5|3cl/2 . (161) 
4 5 

These properties will be extremely useful for checking the 
numerical accuracy of the numerical calculations. 

7 CONCLUSIONS 

We studied the Sunyaev-Zeldovich effect for the clusters of 
galaxies in the Thomson approximation. In Section 2, we 
started with the Lorentz invariant expression of the Boltz- 
mann equation for the electron-photon scattering in the 
CMB frame. Assuming that the initial CMB photon dis- 
tributions are spatially homogeneous both in the CMB and 
CG frames, we derived the formal relation for the rate equa- 
tions between the C MB and CG frames, which agreed with 
IChluba et all (|2012l ). We extended the formula to the rela- 
tion between the CG and OBS frames. The formal relations 
for the rate equations among three Lorentz frames have been 
established under the assumption of the spatial homogeneity 
for the initial CMB photon distributions. 

In Section 3, we investigated the rate equation for the 
photon distribution function in the CG frame. We intro- 
duced three redistribution functions Pe, c {s c ) (£ = 0,1,2) 
and expressed the rate equation in terms of these functions. 
We also reexpressed the rate equation in the operator rep- 
resentation, and the explicit forms were calculated in the 
Fokker-Planck expansion approximation. 



In Section 4, we transformed the rate equation in the 
CG frame to the CMB frame with the formula established 
in Section 2. We co mpared the present result with the pre- 
vious works such aslNozawa. Itoh and Kohv ama (1998j) and 
IChlubaet aI1(|2012l ). We found that the present work agreed 
with these works co mpletely when we multiplied the factor 
1/jc to the result of lNozawa. Itoh and Kohvamal (1 19981 ). 

In Section 5, we transformed the rate equation in the 
CG frame to the OBS frame. The obtained expression had 
advantages compared with previous works. First, the photon 
number conservation was manifestly realized in the present 
form. Secondly, in the limit of /3 a — > /3 C the present expres- 
sion was correctly reduced to the form in the CG frame as 
it should. 

Finally, we derived analytic expressions for the redistri- 
bution functions Pi , c (s c , /3) for I = 0, 1, 2 in Section 6. With 
these formulas, the 5-dimensional integrals for solving the 
rate equations were reduced to the 2-dimensional integrals, 
which will save the numerical computation time drastically. 
These analytic expressions will be extremely useful for the 
analysis of the forthcoming observations. 
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APPENDIX A: ELECTRON NUMBER 
DENSITIES IN DIFFERENT LORENTZ 
FRAMES 

In the present appendix, we discuss relations of the electron 
number densities (JV°, N e and N°) among three different 
Lorentz frames. Let us first define the isotropic electron dis- 
tribution function f c (E c ) in the CG frame as follows: 

2 lw? f ^ )=N - (A1) 

where (E c ,p c ) is the electron four momenta and N£ is elec- 
tron the number density. The electron distribution p etC (E c ) 
used in the main text is related by f c (E c ) = N^Tv 2 p e ,c(E c ). 
It is assumed that the CG is moving with a velocity /3 C with 
respect to the CMB frame. Then, we define the electron 
distribution function f(E) in the CMB frame as follows: 

2 /(S^ /(s)=7v - (A2) 

where (E, p) is the electron four momenta and N e is electron 
the number density. Finally, we define the electron distribu- 
tion function f (E ) in the OBS frame as follows: 



d?p 
(2tt) 3 



f (E ) = iV e ° : 



(A3) 



where (E ,p ) is the electron four momenta and N° is elec- 
tron the number density. It is assumed that the OBS is mov- 
ing with a velocity j3 with respect to the CMB frame. 

Let us first discuss the relation between N e and N£. 
The electron energies E and E c are related by the Lorentz 
transformation: 



E = 7e 1 + 



Pc • Pc 



E c 



(A4) 



The electron distribution function is Lorentz invariant, 
namely, 



/(£) = f c (E c ) . 

Inserting Eq. (|A5[) into Eq. (|A2[) . one has 

d 3 p 



(A5) 



N, 



(2n} 
d 3 p c 



(2*0 



t7c | 1 + ^ ] f e [E c 



E c 



(A6) 



The last equality in Eq. (|A6[) was derived by the following 
reason. The Lorentz invariant phase space volume element 
is expressed by 

therefore d 3 p/{2ix) z receives the same Lorentz transforma- 
tion as that of E, namely, Eq. ()A4[) . The second term in 
Eq. (IA6[) becomes zero by the integral over the solid angles, 
because f c (E c ) is assumed to be isotropic. Finally, one has 



(A8) 



The electron number density in the CG frame is Lorentz 
contracted compared with that in the CMB frame. 

We now discuss the relation between N° and N£. The 
electron energies (E and E) and momenta (d 3 p and d 3 p) 
are related by the Lorentz transformations: 



l: - -. ! i - ^ l • 



(A9) 



(A10) 



where 7 = l/Wl — /3%. In order to Lorentz boost the right 
hand side of Eq. (|A10jl into the CG frame, we first decom- 
pose p into 



P = P\\Pc + P±P± , 



(All) 



where p« = j3 c ■ p and p±fi± = p — (j3 c ■ P)$c- The vectors /3 C 
and /3j_ satisfy the following conditions: 



Then, one obtains the following transformations 

Pc 'P. 



P\)Pc _ IcPc ■ {pc + PcE c ) p c 



E r 



+ Pc)P, 



E 



7c | 1 + I E, 



1 + 



P±P± 
E 



Pc±Pl 



pc_ 

E c 



Pc ■ Pc 

E c 

Pc ■ PcPc 

E c 



7c U + 



Pc ■ Pc 

E c 



(A12) 
(A13) 



(A14) 



(A15) 



Thus, one has 

Pc - Pc 



p_ 

E 



E, 



+ Pc)P 



Pc _ Pc ■ PcPc 

E c 



E c 



1 + 



Pc ■ Pc 

E c 



1c 1 + 



Pc ■ Pc 

E c 



(A16) 



Note that Eq. (|A16|I is correctly reduced to Eq. (fig)) in the 
massless limit. Inserting Eq. (|A16[) and 



d 3 p = 7c(l + ^Jd 3 Pc 

into Eq. (|A10[I . one finally obtains as follows: 



(A17) 



d Po — ~(oJc 



1 - Pa ■ Pc + < Pc ~ Po ■ PcPc 



--{Po-Po-PcPc) 



d 3 p c . (A18) 



One can rewrite Eq. (| A3|) with Eq. (|A18[I and the Lorentz 
invariance of the electron distribution function 

fo(E ) = f c (E c ) (A19) 

as follows: 

K = loic (1 - Po ■ Pc) K , (A20) 

where the terms in the { } brackets in Eq. (|A18|) were 
dropped by integral over the solid angles, because f c (E c ) is 
assumed to be isotropic. Finally, it is clear from Eq. (|A20[I 
that 



AT e =7c(l-/3c)AT e c = iV e c 
in the limit f3 = f3 c . 



(A21) 
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APPENDIX B: NORMALIZATIONS OF THE 
REDISTRIBUTION FUNCTIONS 

Bl Normalization of Po, c (s c ) 

As described in Section 6.1 in detail, we use the ER frame 
to calculate the normalization of Po,c(s c ). It is expressed as 
follows: 

dSc Po,c(s c ) 



dp c plpe,c(E c ) / duo / dy! 



64tt 2 







-1 J -1 

1 



fl'K f It; 

X / d ^J ^?(l-^0)3 



Xa,o ■ 



where 



Xa,o = 1 + cos 6o . 



(Bl) 

(B2) 
(B3) 



The integrals of 0fe c and </> fc / can be done immediately, which 
give 

, s a / #fc c / #&' c (l + cos 2 6o) 
l 27r i Jo Jo 

= 1 + MoMo 2 + 2 ( 2 ~ Mo) (l - Mo 2 ) ■ 

Then, fi' and /io integrals can be also done as follows: 

dfi' 1 1 + filfi' 2 + - (l - /Lto) (l - Mt) 2 ) | = | 



1 

2^ 



(1-/W 
Finally, one obtains 



= 1. 



ds c P ,c(s c ) = / dp c p c p e , c (E c ) = 1 . 
/o 



(B4) 

(B5) 
(B6) 

(B7) 



Finally, we calculate the normalization of the redistri- 
bution function Pnik,i( s )- Inserting D z = into Eq. (|100[) . 
one has 

Onik,i(0) = Co(0) . (B12) 
Thus, the normalization of Pnik,i(s) is given by 



dsPmK,i{s) — / ds Po, c (s) — 1 . 



(B13) 



Equation (|B13[I is consistent with the normaliza- 
tion of the redistribution funct ion Pi,k(s) defined in 
iNozawa. Kohvama and It oh (2009), namely, 



dsP ltK (s) = 1. 



(B14) 



B2 Normalizations of Pe, c {s c ) in the Fokker-Planck 
expansion approximation 

An alternative approach to calculate the normalizations for 
Pe,c{s c ) is to use the operator representation formalism stud- 
ied in Section 3.3. Inserting D z = into Eq. (|50|) . one ob- 
tains the normalizations for Pi >c {s c ) as follows: 



(B8) 



ds c p ttC (s c ) = i + e> £ (o), 

where Oe(0) are given up to 0(6%) by 
0o (0) = 0, 



(B9) 



Oi (0) = 



5 5 140 28 
23363 „ 5 
+ ^48~ 6e ' 



(B10) 



02 (0) 



9 3 183 
10~5 ee + V 
26961 o5 



429 3 2559 4 
w e + —9 e 
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(Bll) 
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